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Abstract. Force and conductance were simultaneously measured during the 
formation of Cu-Cgo and C6o _ C6o contacts using a combined cryogenic scanning 
tunneling and atomic force microscope. The contact geometry was controlled with 
submolecular resolution. The maximal attractive forces measured for the two types 
of junctions were found to differ significantly. We show that the previously reported 
values of the contact conductance correspond to the junction being under maximal 
tensile stress. 
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1. Introduction 

When a molecule is contacted by electrodes to measure the conductance of the molecular 
junction, new bonds are formed and significant forces may arise. These forces affect the 
atomic-scale junction geometry, which is crucial for its transport properties jU El El HJ [5] . 
Current and force can be measured simultaneously using a combination of scanning 
tunneling microscopy (STM) and atomic force microscopy (AFM). Such measurements 
were carried out for metallic contacts (2j EJ El El E]. Related data was reported for 
contacts to single molecules in a liquid environment [101 [TT] and for molecules on a 
metal surface [12]. However, the exact contact geometry was not accessible. 3,4,9,10- 
perylene-tetracarboxylicacid-dianhydride was probed in ultrahigh vacuum using AFM to 
controllably lift the molecule [13] ■ A bimodal distribution of conductances was observed 
and suggested to reflect two distinct bonding geometries. As to controlled molecule- 
molecule contacts, experimental results are few. The conductance of Cqq-Cqo contacts 
was measured by attaching a Cqo molecule to an STM tip and approaching a second 
molecule in a monolayer on Cu(lll) [H]. The force between a metal tip and Cqq 
molecules in double layers on Cu(lll) was addressed with AFM [16l [17]. While close 
distances well into the repulsive range were explored, the corresponding conductances 
[T8] were significantly lower than in the STM work of Ref . [H] . A possible origin for this 
difference may be different geometries of the contact between the tip and the molecule 
in these experiments. Atomically sharp electrodes were shown to act as bottlenecks for 
charge injection into Ceo [151 fT9] . While tips had been intentionally flattened to firmly 
attach a molecule in Ref. [H], the tip used in Ref. [16] presumably was atomically 
sharp. Another possible reason for reduced conductance is foreign material at the tip 
apex. Here, we present low-temperature force and conductance data for the controlled 
formation of Cu-Cgo and Cqq-Cqq contacts. The orientations of the molecules at the tip 
and the surface were determined from STM imaging. The elasticity of both contacts is 
analyzed and compared with density functional theory (DFT) calculations. 

2. Experiment 

Our experiments were performed with a homebuilt STM/AFM in ultrahigh vacuum at 
a temperature of 5 K. Clean Cu(lll) surfaces were prepared by repeated sputtering and 
annealing cycles. Submonolayer amounts of Cgo were then deposited onto the sample 
by sublimation at room temperature. Subsequent annealing to ~ 500 K led to a well- 
ordered 4x4 structure of C 60 [201 Ell [221 123] ■ After additional sublimation of small 
amounts of C 60 onto the cooled sample, isolated C 60 molecules were found on both the 
Ceo islands and the bare Cu substrate [21]. A Ptlr tip was attached to the free prong 
of a quartz tuning fork oscillating with an amplitude of (3 ± 0.2) A at its resonance 
frequency of ~ 28 kHz. The tip was covered with Cu by repeated indenting into the 
substrate until submolecular resolution was achieved. The vertical force F acting on the 
tip at the point of closest tip approach was calculated from the measured frequency shift 
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Af{Az) (shown in SI) as a function of the vertical piezo displacement Az using the 
formalism of Sader and Jarvis [25] . Due to the limited bandwidth of the transimpedance 
amplifier, the current recorded with the oscillating tip is averaged over the entire range 
of oscillation. The non-averaged conductance G(Az) was calculated using the method 
of Sader et al. [26] , which recovers the instantaneous current at the point of closest tip 
approach. The bias voltage V was applied to the sample. Further experimental details 
can be found in the supplementary data. 

We note that the intrinsic energy dissipation of the tuning fork did not change 
significantly during the contact formation. In addition, STM images taken before and 
after the contact measurement showed no changes. These facts suggest that no inelastic 
deformations of tip or molecule occurred. 

3. Cu— Ceo contacts 

Figure [T^a) displays a typical constant-current STM image of a C 60 island used for 
contact measurements with Cu-covered tips. The island comprises two domains which 
differ by an azimuthal rotation of the h-C^o molecules by 60° (hO and h60). Ceo molecules 
adsorbed on Cu(lll) with either a hexagon (h-C^o) or a pentagon (p-C^o) facing up give 
rise to distinctly different patterns in the STM image [27]. Fig. QJb) shows the different 
orientations of the molecules as viewed from the tip position. Double bonds separating 
two hexagons (6:6 bonds) are marked by red bars. Contact data recorded with a Cu tip 
above the center of a h-C^ molecule are shown in Fig. Hfc). Both the total interaction 
force F (solid line) evaluated at the point of closest approach to the sample and the 
instantaneous conductance G (dots) are displayed versus the piezo displacement Az. We 
first focus on the force, which is shown over a wider range of displacements in Fig. CE^d). 
To minimize the electrostatic force, which results from the contact potential difference, 
a bias voltage of V = 0.1 V was applied during the contact measurement [28], [29], [30] . At 
large tip-sample distances, F reflects the long-range van-der-Waals force between the 
tip and the sample. It can be approximated by a power law Fi(Az) = a(Az — Az) b [31] 
with typical fit parameters a = — 5.5nN/A , b = —2.3, and Azq = 7.5 A (fit range: 
Az < A). The fit is shown in Fig. QJd) as a dashed line. The exponent close to 2 
indicates an effective sphere-plane geometry of the junction. 

The short-range force F s (Az) [dashed-dotted line in Fig. Hfd)], which only acts 
on the atoms in the immediate vicinity of the molecular junction, is estimated as 
F s — F — F,. It is attractive for large tip heights, reaches a minimum at Az » 4.8A, 
and finally becomes repulsive for Az > 5.5 A. The fitted van-der-Waals force at contact 
(« 0.5 nN) is consistent with estimates for a sharp tip [32]. We note that the total 
force F and the short-range force F s exhibit maximal attraction at nearly the same 
position Az. In other words, long-range forces do not significantly affect the point of 
maximal attraction. However, they do affect the value of the maximal attraction. Using 
different tips we found it to scatter between 1.5 nN and 2.2 nN. We ascribe the origin 
of these significant short-range forces to the chemical bond formation between tip and 
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Figure 1. (a) Pseudo-three-dimensional illuminated constant-current STM image 
(1.7V, 0.55 nA, 3.8 x 5nm 2 ) of a Cgo island showing three different molecule 
orientations (hO, p, and h60). (b) Sketches of three orientations of Ceo on Cu(lll) 
as viewed from the tip. Red double-bars indicate the bonds separating two hexagons 
(6:6 bonds), (c) Total force F (solid line) and instantaneous conductance G (dots) 
calculated from the simultaneously measured frequency shift and averaged conductance 
(Fig. SI) using Refs. [25J [26] - A voltage of V = 0.1 V was applied to the sample. The 
oscillation amplitude of the tuning fork was A = (3 ± 0.2) A. Dashed lines define the 
point of contact, (d) Total force F (solid line) over a wider range of piezo displacements 
than in (c), fitted long-range force F (dashed line, fit range : Az < OA, sec text for 
details) and short-range force F s = F — F/ (dashed-dotted line). Az = OA corresponds 
to the position defined by the STM set-point I = 0.55 nA at V = 1.5 V. 



molecule. Interestingly, calculated interaction forces for a Si tip and C 60 on Si(100) 
are in a very similar range (1.4. . . 2.0nN) [32]. The F(Az) data of h-C^o and p-C^o 
were very similar except for a shift of w 0.4 A along the abscissa due to the different 
apparent heights of the molecules. This apparent insensitivity to the detailed bonding 
geometry may be attributed to the high reactivity of 6:6 bonds. It causes the Cu atom 
at the tip apex to laterally relax [15] . As a result, a 6:6 bond is most likely contacted 
independent of the orientation of the molecule. The conductance G(Az) in Fig. [T^c) 
shows a typical transition from tunneling at small Az to electrical contact. To define 
the point of contact, the intersection of linear fits in the transition and contact regime 
is used [Fig. 1(c), dashed lines] [5]. The resulting contact conductance G c ~ 0.2 Gq 
(Gq = 2e 2 /h is the conductance quantum) is in agreement with previous experimental 
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results [5J HU dn]. Comparing the conductance data with the simultaneously measured 
force [Fig. 1(c), solid line] we find that the point of contact corresponds to maximal 
attractive force. The same observation is made for C6o~C6o contacts (see below). A 
similar result has been reported from metal-metal point contacts, where a maximal 
attractive force was measured at G ~ Go [2] • Modeling of metallic contacts also suggests 
that the deformation of the junction is maximal at the point of contact formation 
[33| IM] . Recently, two- level fluctuations of the conductance on a /is time scale have 
been reported from Ceo on Cu(100) at the transition from tunneling to contact for a 
metal-C 60 contact [35] . In the present case of a 4 x 4 structure of C 60 on Cu(lll), 
and at the rather low bias voltages used (0.1 V . . . 0.3 V), such fluctuations were not 
observed [36] . 

4. C60-C60 contacts 

By approaching the tip sufficiently close, a single C 60 molecule was attached to the 
tip apex. The orientation of such Ceo tips was determined by 'reverse imaging' on 
small Cu clusters which had been deposited before from the tip onto the bare Cu(lll) 
surface [El [15]. Constant-current images of such a Cu cluster recorded at V = — 2 V 
reveal the second lowest unoccupied orbital (LUMO+1) of the molecule [13] . Compared 
to normal STM images of the LUMO+1 [Fig. [U^a)], they show a mirror image of the 
molecule [37] . 

The relative orientations of the tip and sample Cqq molecules strongly affect the 
conductance of the junction in the tunneling range. Figure [2(a) shows a Ceo island 
imaged using a Cgo-functionalized tip with a hexagon facing towards the surface (h- 
C$o tip). Similar to Fig. Q^a), the island comprises two rotational domains of h-Ceo 
(hO and h60), as well as a few p-C^o molecules. Owing to the different orientations, 
distinctly different patterns are observed with the h-C^o tip for hO, p, and h60 molecules 
[Figs. E]^b-d)]. For example, the center of h-C^o appears either as a maximum (hO) or 
as a minimum (h60) in the STM image. On p-CqqS, threefold symmetry of the h- 
Cqo tip is clearly discernible, which reflects the 5:6 bonds of the molecule at the tip. 
The symmetries of these patterns can be understood from a convolution of the local 
densities of electronic states (LDOS) of the tip and the sample. At V = 1.6 V electrons 
essentially tunnel from the highest occupied molecular orbital (HOMO) of the h-Ceo 
tip to the lowest unoccupied molecular orbital (LUMO) of the molecule at the surface. 
A two-dimensional convolution of these orbitals according to the orientations given in 
Fig. S2 is shown in Figs. [2](e-g) [37]. It reproduces the experimental data with the best 
agreement obtained for the hO pattern. 

Figure Wfa) displays the force (solid line) and the conductance (dots) measured on 
a h-C eo molecule with a h-C eo tip at an applied voltage of V = 0.3 V. Compared to the 
data from a Cu tip [Fig. DD(c)], the maximal attractive force is smaller by a factor of 4. 
In experiments with different C^o tips, this maximal attractive force varied from 0.3 nN 
to 0.4 nN. In part, this scatter may be attributed to the uncertainty of the lateral tip 



Force and conductance during contact formation to a Cqq molecule 



6 




Piezo Displacement Az (A) 



Figure 2. (a) Constant-current STM image (1.6 V, 0.55 nA, 8 x 4nm 2 ) of a Cgo 
island measured with a h-Cgo tip. (b-d) Constant-current closcup images of molecules 
with orientations hO, p, and h60, respectively, (e-g) Simulated images showing the 
convolution of the modeled local density of states (LDOS) of the HOMO of the Ceo 
at the tip and the LUMO of the Ceo at the surface for hO, p, and h60 orientations, 
respectively [37]. (h) Instantaneous conductance G (dots) and total force F (solid line) 
for the h-CsQ-h-Ceo contact calculated from the simultaneously measured frequency 
shift and averaged conductance (Fig. SI) using Refs. [25j[26]. An oscillation amplitude 
A = (2.5 ± 0.2) A was used. A voltage of V = 0.3 V was applied the the sample. The 
dashed line indicates an exponential fit to the conductance for small Az. 

position, which we estimate to be ~ 10% of a C 60 diameter. The conductance measured 
for a C60-C60 contact [Fig. EJ^h), dots] starts to deviate from a purely exponential 
behavior (dashed line) at a piezo displacement close to the position of the maximum of 
the attractive force. The conductance at this point (~ 3 x 10~ 3 Go) is approximately 
two orders of magnitude smaller than with a Cu tip [T3j. When we approached the tip 
further towards the surface until the total force exceeded 0.5 nN, a rotation of the Ceo 
molecule at the tip occurred. 
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Figure 3. (a) Conductance G versus Az during the formation of Cu-/i-C6o (crosses), 
and h-Ceo-h-Ceo (dots) contacts. The solid lines are fits of G^f° rm to the measured 
conductance G using Eq. (P with k = 19.4 Nm" 1 for Cu-/i-C 60 and k = 18.9 Nm" 1 for 
h-Ceo-h-Ceo. More details can be found in appendix B. (b) The junction deformation 
d extracted from the fits to G(Az) for Cu-/i-C6o (solid line) and h-Ceo-h-C^o (dashed 
line). 



5. Comparison of elasticities of Cu— Ceo and Ceo-Ceo contacts 

The forces at the junction cause atomic relaxations which affect the conductance. Below, 
a simple model is used to estimate the deformation of the junction from the measured 
conductance and force data. First, the conductance of a rigid junction is calculated as a 
function of the piezo displacement, G^^(Az) [HJ |37], which is shown in Fig. IB1I The 
junction deformation d is approximated by a linear relation d = F(Az)/k, using the 
experimentally determined force F(Az). We then obtain the theoretical conductance of 
the deformed junction, G t ^° Tm (Az), which depends on the stiffness of the junction k: 

From a fit of G^° Tm (Az) to the measured G(Az) data [Figs.DD[c) and[2]^h)] the stiffness 
k is obtained. Figure |3^a) shows fits (solid lines) for Cu-/i-Ceo (crosses) and h-Cm-h-Cfto 
(dots) contacts. Cu-p-Cgo data (not shown) yield similar results. From measurements 
with different tips we determined elasticities k Prf (16 — 37) Nm -1 for Cu-Ceo and 
k w (13 — 24) Nm" 1 for Cqo-Cqo. The extracted deformation d shown in Fig.^b) 
corresponds to a reduction of the tip-molecule distance of ~ 0.9 A for the Cu-Cgo 
contact (solid line). For the C 60 -C 6 o contact (dashed line), d is smaller (w 0.26 A) and 
the transition from tensile to compressive deformation occurs within the Az range that 
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was accessible in our experiment. While the deformations are smaller than the values 
reported for metal contacts [331 [38], they still significantly affect the conductance. 

The values for k may be interpreted in terms of the elasticities of the components of 
the junctions. DFT calculations taking into consideration several atomic configurations 
were used to estimate elasticities of tip and sample (see appendix A). As summarized 
in Tab. IA1| we find that a metallic Cu tip can be characterized by a stiffness in the 
range k^ ~ (45 — 55) Nm _1 depending on its atomistic structure and on details of the 
calculational scheme. For the sample we find k s ^ plc ~ (H2 — 129) Nm" 1 for h-C^ or 
p-C 60 on reconstructed Cu(lll) [23]. Combining tip and sample elasticities in series, 
we thus estimate k e s ~ (32 — 39) Nm" 1 for a Cu-Ceo contact. Similarly, for a /i-Cgo-tip 
we find fc^o « (43 - 81) Nm" 1 , which leads to k cS « (31 - 50) Nm" 1 for a h-C m -h-C 60 
contact. For both contacts, the experimentally determined elasticity is softer than the 
calculated one by a factor of 2. We attribute this difference to two main factors: First, 
the elasticity calculations can be considered as upper bounds as only a finite number 
of degrees of freedom are taken into account (see appendix A). Second, the elasticity 
estimates above do not take into account the softening of the springs close to contact 
due to the formation of chemical bonds between tip and sample. 

It is instructive to compare the obtained junction stiffness values with that of 
an isolated Cqq molecule. From our DFT calculations we find that squeezing a Cqq 
molecule between two opposite hexagons corresponds to a spring constant of 222 N/m 
(appendix A), i.e., the elasticity of a Cqq is significantly stiffer than the molecular 
junctions considered in this study. The junction deformation therefore mainly involves 
the metal-molecule bonds and the STM tip. 

6. Conclusions 

In summary, simultaneous force and conductance measurements for Cu-Cgo and Cgo^Cgo 
contacts have been performed. Angstrom-scale deformations of the contacts and effective 
stiffness extracted from the experimental data agree with elasticities determined with 
DFT calculations. We find that the maximal attractive force measured at a C6o~Cgo 
contact is 4 times smaller than in a Cu-Cgo junction. Moreover, the force data reveal 
that previously reported contact conductances correspond to geometries in which the 
junctions are under maximal tensile stress. 
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Appendix 

Appendix A. DFT calculations 

In order to estimate the elasticities related to the experimental sample and tip sides we 
considered the generic structures shown in Fig. IA1I Calculations were performed using 
the SIESTA [39] pseudopotential density functional theory method with the PBE-GGA 
exchange-correlation functional [30], a 400 Ry mesh cutoff, and a 2 x 2 Monkhorst- 
Pack [31] fc-mesh. The Fermi surface was treated by a second-order Methfessel-Paxton 
[42] scheme with an electronic temperature of 300 K. The basis set consisted of default 
double-zeta plus polarization (DZP) orbitals for C and Cu atoms generated with an 
energy shift of 0.01 Ry. The force tolerance for the structural optimizations was 0.02 
eV/A. The total energies from SIESTA were not corrected for basis set superposition 
errors. Two different lattice constants for the Cu crystal (a = 3.62 A and a = 3.70 A) 
were considered to confirm that the results do not depend sensitively on this parameter. 




Figure Al. Structures considered in the DFT calculations (top and side views): (a) 
Cu adatom on flat Cu(lll), (b) h-C^o on flat Cu(lll), (c) pyramidal Cu tip, (d) h- 
Ceo on reconstructed Cu(lll), (e) /i-Cgo-tip on Cu pyramid, (f) p-Ceo on reconstructed 
Cu(lll). The boxed atoms in the (4x4)-Cu(lll) slab are kept fixed at bulk coordinates 
while the other degrees of freedom are relaxed. The elasticities are estimated from 
the energy increase associated with displacing the top-most atoms (box with arrow) 
perpendicular to the surface film while allowing the remaining degrees of freedom to 
relax. 

The computational procedure consisted of the following steps: (1) Relaxation of 
initial geometry with the boxed atoms in the Cu slab (Fig. IA1I) kept fixed at bulk 
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Figure A2. Energy differences versus displacement z of the topmost atoms for the 
systems shown in Fig. IA1I The relaxed geometry corresponds to z = A whereas 
tensile strain is exerted on the system for z > A. The solid lines are quadratic fits to 
the calculated data points. The corresponding fitted spring constants are reported in 
Tab. El 

coordinates while the other degrees of freedom fully relaxed. (2) Displacement of the 
topmost atoms (box with arrow in Fig. IA1|) perpendicular to the surface film while 
relaxing all remaining degrees of freedom. 

The energy costs associated with the deformations with respect to the displacement 
z of the topmost atoms are shown in Fig. IA2I with z = A being the relaxed junction. 
For z > A, tensile strain is exerted on the system. Quadratic (as well as fifth order) 
fits to these energy differences yield the effective spring constants in Tab. IA1I It should 
be noted that the elasticities in Tab. IA1 1 represent in fact upper bounds because of the 
limited size of the unit cell. In the real system more atoms respond to the pull on the 
topmost atoms, which leads to a smaller effective spring constant. 

To test the accuracy of the SIESTA calculations, selected checks with the VASP 
code based on a plane-wave basis and the projector- augmented wave method (PAW) 
[1311111115] were also performed. We used PBE-GGA gjD], at least 400 eV energy cutoff, 
a 2 x 2 (or 3x3) Monkhorst-Pack jH] fc-mesh, first-order Methfessel-Paxton jl2] scheme 





Figure A3, (a) Energy differences versus hexagon- hexagon distance z for different 
deformations of an isolated Ceo molecule. The solid lines are quadratic fits to the 
calculated data points, (b) Fixing the distance z between opposite hexagons while 
relaxing all other degrees of freedom (black circles), (c) H ff (l) vibrational mode of Ceo 
(compression/expansion, haj = 31 meV, red diamonds), (d) A ff (l) vibrational mode of 
Ceo (isotropic deformation, fku = 60 meV, blue squares). 
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Table Al. Calculated spring constants Fiesta an d &vasp for the systems shown in 
Fig. IA1I using the SIESTA and VASP codes, respectively. The values in parentheses 
are derived from the second derivative of 5th order polynomial fits evaluated at the 
energy minimum. Two different lattice constants (a = 3.62 A and a = 3.70 A) were 
considered for the Cu crystal. For comparison the stiffness of an isolated Ceo from 
Fig. [A3] is included. 



with 0.05 eV smearing width, and 0.02 eV/A force tolerance. As seen in Fig. IA2I and 
Tab. IA11 the two codes yield similar estimates of the elasticities. 

To determine effective elasticities for the combined elasticity of tip and sample the 
springs from Tab. lAll are added in series. In this way, for a Cu-adatom [Fig. IAl( a)] or for 
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a sharp pyramidal Cu-tip [Fig. lAlf c)] in contact with a H-Cqq on reconstructed Cu(lll) 
[Fig. EUd)], we obtain k cS = (l/k tip + 1/k^y 1 « (32 - 39) N/m. In the case of the 
contact between a Cu tip [Fig. lAl( a) or (c)] and a h-Ceo on the flat Cu(lll) [Fig. IAl( b)]. 
an effective elasticity of k cf[ « (29 — 33) N/m is obtained. Due to the reduced number 
of bonds of Ii-Cqq on flat Cu(lll) in comparison with h-Ceo on reconstructed Cu(lll), 
k c s is smaller. Similarly, for the C 60 -tip [Fig. lAlf b) or (e)] in contact with a h-C G0 on 
reconstructed Cu(lll), we estimate k c s ~ (31 — 50) N/m. 

We also analyzed the stiffness of an isolated Ceo molecule with SIESTA, cf. Fig. [A3] 
By controlling the distance between two opposing hexagons while relaxing all other 
degrees of freedom we obtain an effective stiffness of the molecule of k c g = 222 N/m. 
A comparable (but slightly larger) stiffness of fcngi = 306 N/m is obtained when 
considering only a deformation along the characteristic elongation/compression H s (l) 
vibrational mode [Fig. IA3( c)]. Deforming along the isotropic A s (l) vibrational mode 
[Fig. IA3( d)] yields a much larger stiffness (&Agi = 2942 N/m). 

Appendix B. Elasticity model 

The influence of the junction deformation on the conductance is described as follows. 
First, the conductance of a rigid system as a function of the piezo displacement, 
G^^(Az), is calculated for Cu-/i-C 60 and C 60 -C 6 o contacts (Figs. IBl( a) and (b)). 
Az = corresponds to a distance between the topmost tip atom and the topmost 
C 60 -hexagon of 8 A in the Cu-/i-C 60 contact and to a C 60 -C 6 o center distance of 14.9 A 
in the Ceo~C6o contact. Further calculational details can be found in Ref. |14j . 

Next, we assume that the deformation d depends linearly on the experimentally 
determined force, i. e. d — F(Az)/k with k describing the effective stiffness of the 
junction. The theoretical conductance of the deformed junction, G t ^° Tm (Az), depending 
on k is given by: 

Az) = G^S (A* - • (B.l) 

The interpolated calculated conductance G^° im is then fitted to the experimental 
conductances G(Az). As the absolute tip-sample distance was unknown in the 
experiments, an arbitrary shift on the abscissa for G^° d (Az) is allowed. This shift 
is already comprised in Fig. \B1\ so that the Az axes in Fig. IBl( a) and (b) correspond 
to the axis in Fig. 3. 

Fits for the Cu-/i-Ceo and C6o~C6o contacts shown in Fig. 3(a) and to similar data 
calculated for other microscopic tips yield k « (16 - 37) Nrn -1 and k ~ (13 — 24) Nrn -1 , 
respectively. Compared to the values from the DFT calculations in appendix A, 
the effective spring constants qualitatively agree, but are smaller by a factor of 2. 
This deviation is not unexpected since the calculated elasticities are upper bounds. 
Furthermore, the model neglects the formation of chemical bonds between tip and 
sample, which are expected to weaken the effective spring constant. 
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Figure Bl. Interpolated calculated conductance G*]^ as a function of the piezo 
displacement Az for the (a) Cu-/i-Ceo and (b) Ceo-Ceo contact. For the Cu-h-Ceo 
contact, Az = A corresponds to a distance between the topmost tip atom and the 
topmost Ceo hexagon of 8 A. For the Cgo-Cgo contact, Az = A corresponds to a 
C6o~C6o center distance of 14.9 A. 
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1. Experimental details 

Frequency shift-distance curves were recorded as follows. The current feedback was 
disabled and the tip was retracted from the surface by several angstroms. The frequency 
shift Af and the time-averaged tunneling current J avg were recorded while approaching 
the tip towards the surface and during the subsequent retraction to its starting position. 
Curves were acquired with 2000 data points in 60 to 120 s. Since no significant 
differences were observed between the approach and retraction data, modifications of 
the tip apex or the molecular orientation during measurements can be excluded. In 
addition, STM images taken before and after current-distance measurements showed no 
changes. Furthermore, the intrinsic energy dissipation of the tuning fork did not change 
significantly during approach and retraction of the tip. As the cantilever, one prong of a 
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Figure SI. Frequency shift df and time-averaged conductance G avg versus the piezo 
displacement Az for the Cu-Ceo and C6o _ C6o contacts shown in Fig. 1(c) and Fig. 2(h), 
respectively. The oscillation amplitudes are A = (3 ± 0.2) A for the Cu-Ceo contact 
and A = (2.5 ± 0.2) A for the C 60 -C 6 o contact. 

tuning fork was used (spring constant k ~ 1800 N/m) while the second prong was fixed 
on the piezo scanner. Measurements were performed at oscillation amplitudes between 
A = 1 A and A = 3 A. The amplitude A corresponds to half of the peak-to-peak distance 
of the oscillation. The uncertainty of the calibration of the used amplitude is ~ 7%. This 
margin is essentially due to the uncertainty of the calibration of the z-piezo, where ~ 5% 
are routinely achieved for small displacements. The bandwidth of the transimpedance 
amplifier for the tunneling current varies between 7 and 200 kHz depending on the used 
amplifier gain. The current data was further filtered by a numerical low-pass filter with 
a cut-off frequency of 8 kHz. 

1.1. Extraction of the short-range force 

The short-range force F s can be obtained by subtracting the long-range van-der- 
Waals force Fi from the total force F. We use a power law to approximate Fi by 
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Fi(Az) = o(Azq — Az) b . The fit parameters depend on the range of the fit. E.g., for 
the fit range Az < OA we find a = — 5.5nN/A , b = —2.3, and Azq = 7.5 A. For 
Az < -lA and Az < 1 A we obtain a = -6.5nN/A 6 , b = -2.34, Az = 7.8 A and 
a = — 2nN/A , b = —1.9, Azq = 6.3 A, respectively. Figure IS2l shows the influence 
on F s for the different fit ranges. The absolute value of F s can change by a factor of 
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Figure S2. Total force F (dots) and short-range forces F s (solid, dashed, and dashed- 
dottcd lines) versus the piezo displacement Az. The three F s curves are extracted by 
three long-range force fits with different fit ranges. 

2 depending on the fit range. However, the displacement Az where contact occurs is 
only weakly influenced by the different long-range force fits. It shifts by ~ 0.12 A at 
most. Therefore we did not use F s for a quantitative comparison. Rather we compare 
the more robust point of maximal attraction with the point of contact. 

2. Simulation of images recorded with C 60 -tips 

The highest occupied molecular orbitals (HOMO) of Cqq exhibit a large local density 
of states (LDOS) at the 6:6 bonds, which is known from DFT calculations. Motivated 
by this result, we modeled the HOMO of a fo-Cgo-tip by a threefold-symmetric pattern 
of two-dimensional Gaussians [Fig. IS3T a)-(c)]. For simplicity it was assumed that only 
the LDOS of the hexagon which is closest to the surface contributes to the tunneling 
current. The same approach was used to mimic the lowest unoccupied molecular orbital 
(LUMO) of Ii-Cqo and p-C^o on the substrate, which are most intense at the 5:6 bonds 
[Fig. IS37 d-f)]. As images of adsorbed h-C 60 recorded with a Cu-tip reveal the positions 
of the bonds, the relative orientation of Cqq at the tip and the surface were determined 
from the experiment. For p-C^o the situation is less favorable, because the molecule 
appears circularly symmetric when using a Cu tip. A two-dimensional convolution of 
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Figure S3, (a-c) Model LDOS of the HOMO of a ft.-C6o-tip hexagon as viewed from 
the tip. (d-f) Model LDOS of the LUMO of h-C 60 and p-C 60 on the surface, (g-i) 
Two-dimensional convolution of tip and sample LDOS. White hexagons and a black 
pentagon indicate the molecular orientation as viewed from the tip. A, B, and C refer 
to the surface regions defined in Figs. 1(b) and Fig. 2(b). (j) Constant-current STM 
image (— 2 V, 0.26 nA, 1.7 x 1.6 nm 2 ) of a Cu cluster with a Cgo tip. 

these patterns [Figs. 2(e-g) and Figs. I53lf g-i)] reproduces the experimental images of 
h-C^o qualitatively very well. The orientation of the Ceo molecule attached to the tip 
was determined by scanning a Cu cluster, consisting of a few Cu atoms, with the C 60 
tip. Fig. fS3tj ) shows the experimental data (V = — 2 V). This image is rotated by 
180°(about an axis perpendicular to the image plane) compared to that of a Ceo on the 
sample surface scanned with a metal tip. 



